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Abstract 

Meromorphic solutions of autonomous nonlinear ordinary differen- 
tial equations are studied. An algorithm for constructing meromorphic 
solutions in explicit form is presented. General expressions for mero- 
morphic solutions (including rational, periodic, elliptic) are found for 
a wide class of autonomous nonlinear ordinary differential equations. 

1 Introduction 

Autonomous nonlinear partial and ordinary differential equations frequently 
arise in mathematical models describing processes and phenomena in physics, 
chemistry, biology, etc. The problem of constructing exact solutions to au- 
tonomous nonlinear differential equations is intensively studied [THTU]. A 
number of methods and algorithms has been introduced. As a rule these 
methods deal with a priori fixed families of meromorphic exact solutions. 
Consequently, solutions outside given families may be missed. Very often 
various methods yield the same families of exact solutions written in a differ- 
ent way. And as a result many "new" exact solutions appear. [TTHTl] Thus 
an important problem is to classify families of exact solutions. In this paper 
we study the general form of meromorphic solutions to autonomous nonlinear 
ordinary differential equations. The main aim is to present our algorithm, 
which can be applied to construct exact meromorphic solutions in explicit 
form. For a wide class of nonlinear ordinary differential equations we make 
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a conclusion that exact solutions given in our paper are the only possible 
meromorphic solutions. 

This paper is organized as follows. In section 2 we describe our method 
and present explicit expressions for meromorphic solutions. In sections 3 we 
give an example and construct meromorphic solutions of a certain second 
order nonlinear differential equation. 

2 Method applied 

Let us look for nonconstant meromorphic exact solutions of an autonomous 
nonlinear ordinary differential equation 

E[w{z)\ = 0. (1) 

In this expression i£[it?(z)] is a polynomial in w(z) and its derivatives. For 
every solution w(z) of equation ([1]) there exists a family of solutions w(z — 
z ). Without loss of generality, we omit arbitrary constant z . Suppose the 
following. 

Condition I. For solutions of equation ([T]) there are only iV different 
asymptotic expansions corresponding to the Laurent series in a neighbor- 
hood of z = 

Pi (i) oo 

^HE^F + E^ 0<M<*, i = l,...,N. (2) 

k=l k=0 

In this expression pi > is an order of the pole z — 0. If AT = 1 we shall 
omit the upper index. 

Condition II. Substituting w(z) = XW(z) into equation ([rj yields expres- 
sion with only one term of the highest degree in respect of A. 

Note that if w(z) has more than A" poles, then w(z) is periodic. Indeed, 
there exist poles z — z\, z — zi of w(z) such that the functions w(z + zi), 
w(z + z<i) are meromorphic solutions of equation with a pole at z — 0. 
Thus we have the equality w(z) = w(z + z 2 — Zi). 

Meromorphic solutions of equation ([T]) are classified in the following the- 
orems. 

Theorem 1. All meromorphic solutions of equation (JTJ satisfying con- 
dition I with N = 1 and condition II are of the form: 
1) elliptic solutions with the periods 2u±, 2uj 2 

( P (_l)fc c _ fc fjk-2 "j 

w(z) = |E fopi tPteui'Vd + ho, (3) 
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Necessary condition for elliptic solutions to exist is c_i = 0. 
2) periodic solutions with the period T 



^ = l\f}^^^UC4)^- (4) 



T {k-l)\ dz^ 1 j \T ) 

3) rational solutions 



rn 



w W = E^x + E c ^ fc ' m ^°- ( 5 ) 



Z 

k=l k=0 



In 1) and 2) ho is a constant. 

Proof. For any meromorphic function f(z) with the poles {a n } Mittag- 
Leffler found the representation 



f(z)=jr{Pn(j^ r ) -g n (z)\+h(z), (6) 

n=l ^ n ' ^ 

where {P n [(,z — cirj^ 1 ]} are principal parts of Laurent expansions for the 
function f(z) around poles {a n }, {g n (z)} are polynomials, and h(z) is an 
entire function. The polynomial g n (z) is added to the n th term of the sum 
(J6j) in order to provide convergence of the series. A series of the form ([6]) 
is said to be convergent in M C C if only finite number of terms possesses 
poles in M and without these terms the series converges [T5j. Suppose w(z) 
has two or more poles. Then from condition I follows that w(z) is doubly 
periodic (elliptic) or simply periodic. For the elliptic function w(z) with 
Laurent expansion (J2J) expression ([6]) reads 



C-2 

w[z' 

z 

(n,m)^(0, 0) L 



C-2 C_ 2 



— 2nu;i — 2mcj2) 2 (na;i + ma^) 2 



(7) 



\ ^ \ ^ c-fc , 

^ ^ (z - 2nu l - 2mu 2 ) k °' 

k=3 n, rn v ' 



Solution (j2j) can be rewritten in terms of Weierstrass elliptic function p sat- 
isfying the equation 

(p z ) 2 = V - g 2 p - (8) 

Thus we obtain representation for the elliptic solution w(z) in the form ([3]). 
Necessary condition c_i = follows from the theorem for total sum of the 
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residues of an elliptic function in the parallelogram of periods. Writing ex- 
pression (El) for the solution w(z) with period T yields 



w(z) 



-+E 



C-l 



+ 



C-l 



[z — nT) nT 



EE 

k=2 n 



C-k 



- nT)' 



+ (9) 



In the paper [16] meromorphic traveling wave solutions of the Kuramoto - 
Sivashinsky equation were studied. A powerful method based on Nevanlinna 
theory [T7] was developed. From the results of the paper [TH] it follows that 
periodic solutions of equation ([T]) satisfying condition II are rational in e rz . 
Making change of variables s = e rz , y(s) = w(z) in equation ([1]) and then 
substituting y = s n , n > as s tends to infinity and y = s~ l , I > as s tends 
to zero into the result, we wee that an entire periodic function h(z) in (Q is 
a constant since condition II is valid. Further, making use of the formula 



i + E 



1 1 7T ( 7lZ \ 

[(T^T) + nrj = T COt \TJ 



(10) 



we get 01]) from (E]) with h(z) = h . Now suppose w(z) has only one pole, 
then 



w(z 



> = E^ + M 



k=l 



where an entire function h(z) does not have an essential singularity at infinity 
[T6] . This completes the proof. 

Remark. For fixed values of parameters in equation ([1]), if any, there may 
exist only one meromorphic solution (rational, simply periodic or elliptic) 
with a pole at z = 0. 

In the next theorem we generalize results of theorem 1 to the case N > 1. 
Let us call poles of meromorphic solution w(z) with the Laurent expansion 



w 



0/ 



z) (see ([2])) as poles of type i. 



Theorem 2. All meromorphic solutions of equation ([T]) satisfying con- 
dition I and condition II are of the form: 
1) elliptic solutions with the periods 2ui, 2ui 2 



-l) k c ( H d k ~ 2 



W(Z 



(k-l)\ dz k ~' 
Pi 




p z (z) + B L 



4 |_ p{z) - Ai 



EE 

iel k=2 
\- C l l\(p z (z) + Bj) 

+ ^ 2(p(z)-A i ) + \^ 2 (k-l)\ dz^ 



(12) 



p{z) + ho, 
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where p(z) d = p(z; u 2 ), Bf = AAf — g 2 Ai — g 3 . Necessary condition for 
elliptic solutions (fT2|) to exist is 

J2c { l\ + c { ^ =0. (13) 
2) periodic solutions with the period T 



T \hh ^-^ dzk ' 1 \ A-fcot(f) 

vr f A (-l)*- l c<5 d fe - x 1 /7TZN , 

+ r ^-n. ^-i cot It") + h °- 



(14) 



Tl^ (k-l)\ dz^ 1 \ \T 
3) rational solutions 

Pip (io 

Eif+EE77^i+E^ ™so. (i5) 



Pl o (io) Pi „(i) 



fc=l is / fc=l V ' k=0 



In 1), 2) and 3) {A}, {-£>«}, {a,}, {/ij} are constants and / = or / C 
{l,2,...,AT}\{i }, l<i <AT. 

Proof. Let u'(z) be a meromorphic solution of equation ([[]) with poles of 
types j, j G J, JC {1,2,..., iV}. Without loss of generality, we may select 
any i G J and suppose that the point z = is a pole of type io for 10(2:). 
Let us set I = J \ {io}- For N = 1 we see that 1 = 0. Let us consider the 
case N > 1. If there are more than |/| + 1 poles of different types, then w(z) 
is periodic (simply periodic or elliptic). By analogy with the case N = 1 we 
get the expression for elliptic solutions 



w[z) 



A {-i) k c% d k ~ 2 

^ (k-l)\ dz k ~ 2 

i€lU{i } I ieIU{i } k=2 y ' 



(16) 



where ((z) is the Weierstrass ^-function and aj may be taken as zero. Note 
that dj 7^ dj, i, j G / U {io}, i 7^ j are poles of elliptic solutions ( !T6|) in the 
parallelogram of periods built on 2u>i, 2u 2 . The function ((z) has a simple 
pole at the point z = with residue 1. The theorem for total sum of the 
residues of an elliptic function in the parallelogram of periods yields condition 
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ffT5j) . Using the addition formulae [IB] 



CO -a;) = CO) -CK) + o 



p(z - a,) = -pO) - p(a») + £ 



2 p(z) - p { Oi) ' ^ 



def def 

and notation Aj = p(aj), 5j = p 2 (aj), we get equality (1T21) . Further, from 
equation ([8]) we have = — <?2^4i — #3- Note that for the parameters ho 
and ho the following correlation holds 



ho = h - J^ c -iC( a i) - J^c^pOi). 



Similarly to the case N = 1 simply periodic solutions can be written as 



w\z) 



E E (fc-i)! djgfc -i ^ cot — y — J + fao> (is) 

ielU{i } k=l v ' 7 



where again we may set aj = and dj ^ a,j, i, j G / U {«o}j i ^ j are 
poles of simply periodic solutions f jT6|) in a stripe built on T. With the help 
of addition formulae for cot [71*0 — <2j)T -1 ] we obtain expression f|T4|) . where 

A; 7r cot[n diT~ l ]T~ l . For rational solutions we have the equality 



(*) 



E Ett^+E^ m ^ ' ( 20 ) 



— o, . 



where dj ^ dj, i, j G / U {«o}, « 7^ j- Without loss of generality, we set 



a io = in (|20]) to obtain (fT5|) . 

Remark 1. Again for fixed values of parameters in equation ([I]), if any, 
there may exist only one meromorphic solution (rational, simply periodic or 
elliptic) with a pole at z = of type i, 1 < i < N . 

Remark 2. While making classification of meromorphic solutions, we 
should consider all the variants for JC {1,2,..., iV}. 

Now let us discuss the problem of finding exact meromorphic solutions in 
explicit form. The main idea is to expand exact solutions given in theorems 
1 and 2 around their poles and to compare coefficients of these expansions 
with coefficients of series 02). Our algorithm is the following. 

At the first step one constructs asymptotic expansions corresponding to 
Laurent series in a neighborhood of a pole z = 0. 
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At the second step one selects an expression for a meromorphic solution 
(see theorem 1). 

At the third step one expands the meromorphic solution in a neighborhood 
of poles z = 0, z = cii, i G /. 

At the fourth step one compares coefficients of the series found at the first 
and the third steps and solves an algebraic system for the parameters of ex- 
act meromorphic solution. Along with this, correlations for the parameters 
of equation may arise. If the meromorphic solution possesses poles of 
different types, it is convenient to consider all possible Laurent expansions. 
Optimal number of equations in algebraic system equals to the number of pa- 
rameters of the meromorphic solution and equation plus 1. If this system 
is inconsistent, then equation ([T]) does not possess meromorphic solutions 
with supposed expression. 

At the fifth step one verifies the meromorphic solution obtained at the 
previous step, substituting this solution and correlations on the parameters 
of equation (JTJ) into the latter. 

In order to find the Laurent expansions for elliptic solutions (theorem 1) 
in a neighborhood of poles z = 0, z = i G / one takes formula (TTtjj) with 

def def 

a io = and then introduces notation A t = jp(cij), Bi = p z (ai), % G /. The 
parameters of elliptic solutions f[T2l that one needs to calculate are g 2 , #3, 
h , Ai, i G /. 

As far as simply periodic solutions (theorem 1) are concerned, one takes 
formula (TT9j) with aj = 0, constructs Laurent expansions in a neighborhood 

of poles z = 0, z = ai, i G /, and introduces notation L d = n 2 T~ 2 , Ai d = 
7r cot[7rajT _1 ]T _1 , i G I. The parameters of simply periodic solutions ffT4l) to 
be found are the following L, h , Ai, i G /. 

For rational solutions (theorem 2) the parameters to be calculated are 
h , ai, i G I . In order to find the highest exponent m (see ffTB^) ) one may 
construct the Laurent expansion in a neighborhood of infinity. 

To conclude this section we would like to emphasize that our method 
is also applicable to algebraic autonomous nonlinear ordinary differential 
equations with arbitrary coefficients in asymptotic expansions corresponding 
to the Laurent series in a neighborhood of poles. Arbitrary coefficients should 
be added to the list of parameters. However, these equations may possess 
meromorphic solutions of more complicated structure. 
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3 Explicit construction of meromorphic solu- 
tions 

As an example let us find exact meromorphic solutions of the following second 
order autonomous nonlinear differential equation 

w zz + ww z + aw 3 + 5w z + (3w 2 + aw + 7 = 0, cr^O. (21) 

Without loss of generality, we may set 5 = 0, a = —6, and study the equation 

w zz + ww z — 6w 3 + (3w 2 + aw + 7 = 0. (22) 

This equation admits two different asymptotic expansions corresponding to 
Laurent series in a neighborhood of poles. They are the following 

w( ^ = Tz + To + { W + a ) 9 + (375 + 15 + V T + o(N } (23) 



(2)/ x 2,2/3 [20(3 2 , \ z , 3 / 2/3 d 



(24) 

We see that equation ([221 satisfies conditions II and I with N = 2. First let 
us construct simply periodic solutions possessing poles of single type (with 
residue = 1/2 and then with residue cl\ = —2/3). In this case we should 
take exact meromorphic solution in the form 

w(z)=c. 1 VL cot (yZz)+ho, L= f ^. (25) 

Note that we omit arbitrary constant z . Finding the Laurent expansion of 
function ([25]) at the point z = 

„(,) = f=i + ho - <±L Z - ^ _ + (26) 



and comparing coefficients of this series with coefficients of expansion ([23]) 
(and then fl2^|) ). we obtain parameters L, ho and correlations on the parame- 
ters a, (3, 7. Substituting obtained meromorphic solutions into equation f[2"2"j) . 
we make sure that this equation indeed possesses meromorphic solutions of 
the form ([25]) . For the Laurent series (1231) these solutions can be written as 

4)4cot(^UA T = lo(-—* V /2 7r(27) 

v ; 2T \ T J 20 V H/3 2 + 200ay 
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6 



w(z) 





w(z) 



1 



-4 



(a) a = -10, /3 = 10, z = 0, T = V3vr/3 (b) a = 10, /3 = 10, z = 0, T = V93vri/31 

Figure 1: Family of solutions (127|) . 



provided that the following correlation holds 

(/3 2 + 25a)/3 
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375 



(28) 



For the Laurent series (1241) we get the family of meromorphic solutions 



2tt J n(z-zo) \ 2(5 

w(z) = — cot < — > H , 1 = 22 , 

3T 1 T J 33 V 20(3 2 + 363a 



and conditions on the parameters of equation 
(4/3 2 + 121a)/3 



1/2 



7T (29) 



7 



2662 



(30) 



Note that in expressions (127|) . (|29j) the period T depends on the parameters 
a and /3 and can be real or complex. For the functions with purely imaginary 
period T = ir, r e IR the following correlation holds 



-cot ( Y j = 7 C0t HvJ- 



(31) 



Now let us construct exact meromorphic solutions with poles of different 
types. The general expression for this solution is the following (see expression 
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-2 



w(z) 



-2 
-2 - 

-4- 



(a) a = -10, p = 10, z = 0, T = 22tt/V815 (b) a = 10, /3 = 10, z = 0, T = 227ri/\/2815 

Figure 2: Family of solutions (|29|) . 

dH]) with i = 1, I = {2}, a 2 = a) 

wO) = cl^VIcot jv 7 !^} + cl 2 jVZcot {vT(z - a)} + h . (32) 

Again we omit arbitrary constant z . Without loss of generality, the point 
z = is a pole of the first type, while the point z = a is a pole of the second 
type. Constructing the Laurent series for solution (132]) in a neighborhood of 
the points z = 0, z = a, we obtain 



.(1) / Jl) r \ 

^ + / i o-A C L1-[cL 2 J(L + A) + ^j 



z + o(|z|), 



(2) / (2) r 

_^zL_ + ho + A c « - cS(£ + A) + ^- I (2 - a) + o(|* - a|), 
2 — a \ 3 



(33) 

where we set A d = y/Lcot y/~La. Comparing the first of these series with the 
expansion (123 j) and the second — with expansion (j24p , we find the parameters 
ho, T, A and conditions on the parameters a, (3, 7 of equation (l2"2l . Thus 
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we have 



P (Vcot<j 


' 0(z-z o ) 1 

. no J 




61/3 


165 ^cot < 


f /3(js-«o) " 

I 110 . 


} + 7) 


660 



w(z) = £- cot ( ^ 7/ o) }- \ \ - I ', + (34) 



220 [ 110 

The family of meromorphic functions (134"]) solves equation ( 122]) provided that 



94/3 2 31/3 3 



a = -7T7^' 7 



3025 ' 33275 
Expression (jMj) can be rewritten as 



(35) 



„ M = JL cot (gi^l I _ JL cot ( g(£ - gO + 1 , og 1 1 + gig, (36) 

v ; 220 \ 110 J 165 \ 110 2 & 3j 660 V ; 

Meromorphic solutions (I34p with z — possess two poles on real axes, one 
- at the origin, another one - at the point a = — 551og(4/3)//3, /3 G R. 
All other poles lie in the complex domain since the period of solutions is 
T = HOttV/3. 

Now let us construct rational solutions of the equation in question. Sub- 
stituting w = z m as z tens to infinity into equation f[2"2~]) we see that m — 0. 
Note that at certain conditions on the parameters a, /3, 7 expansions fl2"3"|) . 
f[2"4"|) terminate. As a result we find two families of rational solutions and 
correlations for the parameters of equation (122]) 



1 /3 H/3 2 

= — + — , a = ~-^-r, 7 



2(z-z ) 20' 200' 1000 

2 2/3 _20^2 4/3^ 1 j 

^ ~~3(^-^) + 33' "~~1363~' 7 " 3993" 

If we try to find rational solutions with poles of two types, we obtain an 
inconsistent system of algebraic equations. Consequently, rational solutions 
of the form 

(1) (2) 

w(z) = + — ^ + ho. (38) 

z — zq z — a — zq 

do not satisfy equation (122]) . 

Equation (122]) does not posses elliptic solutions, since necessary condition 
(TT3|) does not hold. 

We have studied all the expressions for meromorphic solutions given in 
theorem 2. Hence, we have found all families of nonconstant meromorphic 
solutions of equation (1221 . Plots of these solutions for certain values of the 
parameters are given in figures [1] — [3] 
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4 - 

2 - 


4 - 

\ w(z) 

V 2 " 


\ / 








-4 -2 




A 

-2 - 
-4 


2 4 -4-2 




2 4 

z 


-4 - 





(a) /3 = 10, z = 0, T = llvri (b) /3 = -10, z = 0, T = 11m 



Figure 3: Family of solutions fl34|) . 

4 Conclusion 

In this paper we have presented a method for constructing exact solutions 
of autonomous nonlinear differential equations. Our approach is simple in 
application since one can introduce parameters of exact solutions in such 
a way that for many equations linear correlations on the parameters arise 
at each step. Our method generalizes several other methods (for example, 
exp-function method, tanh-function method [TU] and their extensions and 
modifications [20] )• A great drawback of these methods is an a priori as- 
sumed expression for an exact solution. As a consequences solutions outside 
stated family may be lost. Our approach takes into account the structure 
of singularities and possible Laurent expansions for exact solutions. Along 
with this, we have obtained the general form of meromorphic solutions for a 
wide class of autonomous nonlinear ordinary differential equations. 
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